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Electron spin qubits in silicon, whether in quantum dots or in donor atoms, have long been consid¬ 
ered attractive qubits for the implementation of a quantum computer due to silicon’s “semiconductor 
vacuum” character and its compatibility with the microelectronics industry. While donor electron 
spins in silicon provide extremely long coherence times and access to the nuclear spin via the hy- 
perfine interaction, quantum dots have the complementary advantages of fast electrical operations, 
tunability and scalability. Here we present an approach to a novel hybrid double quantum dot by 
coupling a donor to a lithographically patterned artificial atom. Using gate-based rf reflectometry, 
we probe the charge stability of this double quantum dot system and the variation of quantum 
capacitance at the interdot charge transition. Using microwave spectroscopy, we find a tunnel cou¬ 
pling of 2.7 GHz and characterize the charge dynamics, which reveals a charge T 2 * of 200 ps and a 
relaxation time Ti of 100 ns. Additionally, we demonstrate spin blockade at the inderdot transition, 
opening up the possibility to operate this coupled system as a singlet-triplet qubit or to transfer a 
coherent spin state between the quantum dot and the donor electron and nucleus. 


I. INTRODUCTION 

Interest in donor-based spin-qubits in silicon has been 
motivated by their exceptionally long electron spin co¬ 
herence times, exceeding one second in isotopically en¬ 
riched 28 Si ID- Additionally the donor electron spin can 
be a gateway to access the donor nuclear spin, which has 
longer coherence times [2], even at room temperature j3] 
and the potential to serve as a quantum memory [4]. 
Moreover, the single shot read-out of single electron [5] 
and nuclear [6, spins, a milestone for donor-based quan¬ 
tum computing, has been recently demonstrated in na- 
noelectronic silicon devices. 

On the other hand, artificial atoms such as electro¬ 
statically defined quantum dots offer complementary ad¬ 
vantages as qubits, notably in their tunability [7], flex¬ 
ible coupling geometries [8], and opportunities for fast 
electrical control of spin [9]. It is therefore attractive to 
investigate the possibility of hybrid architectures which 
brings together the advantages of these two systems by 
coupling a quantum dot to a donor atom [T0l [IT] . Such 
a double dot could take advantage of fast spin manipula¬ 
tions using gate voltage to form a hybrid singlet-triplet 
qubit coupled to the long-lived quantum memory offered 
by its nuclear spin. In addition, this hybrid architecture 
could be used to create spin buses with quantum dot 
chains to mediate quantum information stored in donor 
qubits over long distances [12 . 

In this work, we present measurements taken on a 
silicon nanowire transistor indicating a coupled system 
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formed of a single phosphorus atom and a quantum dot. 
We show that this hybrid system behaves as a double 
quantum dot, and characterise the high-frequency ad¬ 
mittance of the coupled system by RF gate-based sens¬ 
ing. Quantum capacitance changes at the interdot charge 
transition (ICT) and microwave (MW) spectroscopy al¬ 
low us to characterize the tunnel coupling and the charge 
dynamics. Finally, we demonstrate spin blockade effects 
between singlet and triplet states by applying a magnetic 
field which changes the resonator response at the ICT. 
These results demonstrate the potential of a donor-dot 
system as a new singlet-triplet qubit. 


II. CHARACTERIZATION OF THE DEVICE 

A. Device fabrication and measurement setup 

The nanowire transistor device, sketched in Fig. [lja), 
is completely CMOS industry compatible and is fabri¬ 
cated from an SOI substrate with a 145 nm buried oxide 
and a silicon layer doped with phosporus at a concentra¬ 
tion of 5x 10 17 cm -3 . The doped silicon layer is etched to 
create a 200 nm long and 30 nm wide nanowire by means 
of deep-UV lithography. A 30 nm wide wrap-around 
top gate is defined using a SiO2(0.8nm)/HfSiON(1.9nm) 
stack for the gate dielectric followed by TiN(5 nm)/poly- 
Si(50 nm) as the top gate material. The source and drain 
are self aligned and formed by ion implantation after the 
deposition of 20 nm long SisN 4 spacers. 

Measurements are performed at the base temperature 
of an Oxford Instrument Triton 200 cryogen-free dilution 
refrigerator (30 mK). High sensitivity charge detection is 
achieved by radiofrequency reflectometry on a tank cir- 
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FIG. 1. (a) A sketch of the device cross-section shows quan¬ 
tum dots localized in the top corners of the nanowire and a 
dopant located deeper in the channel, (b) In our rf reflec- 
tometry measurements, an rf signal is sent through a direc¬ 
tional coupler to a resonator made from the device capac¬ 
itance, parasitic capacitance and surface-mounted inductor. 
The reflected signal is amplified at low temperature, and de¬ 
modulated using a reference signal to give V r f. (c) The deriva¬ 
tive of the demodulated resonator response as a function of 
Vt g and Vb g gives the charge stability diagram. Transitions 
with weak Vb g dependence are attributed to changes in the 
charge occupancies of the two distinct (and uncoupled) cor¬ 
ner dots. The charge transition indicated with a black arrow is 
attributed to a phosphorus atom in the bulk of the nanowire. 


cuit composed of a surface mounted inductance (390 nH), 
a parasitic capacitance to ground (0.75 pF) and the de¬ 
vice capacitance between the transistor top gate and the 
channel [15] . Rf-reflectometry is performed close to res¬ 
onance frequency (294MHz), with the input power set 
to —85 dBm and the reflected signal amplified by a low 
noise cryogenic amplifier anchored at 4K. The signal is 
further amplified and demodulated at room temperature 
as shown in Fig. [ljb). An on-board bias tee is used to 
apply both DC and RF voltages on the top gate. The 
undoped silicon substrate is activated by flashing a sur¬ 
face mounted blue LED to generate free carriers and can 
then be used as a back gate m- 


B. Charge stability diagram 

We first characterize the device charge stability dia¬ 
gram as a function of top gate {Vtg) and back gate volt¬ 
ages (Vbg). Rf-reflectometry is used to detect a change in 
the sample impedance due to a dissipative or dispersive 
event occurring below the top gate. As a consequence 
it can probe a change of resistance, due for instance to 
a charge tunneling between the source and a localized 
state [HI US], or a change of quantum capacitance in¬ 
duced by a charge tunneling between two quantum dots 
Dana. Fi g . c) presents the filtered numerical deriva¬ 
tive of the reflected resonator signal amplitude as a func¬ 
tion of V tg and Vb g . It shows a set of charge transitions 
with a small dependence on Vb g corresponding to single 
charge tunneling from the source or drain to quantum 
dots localized below the top gate. We attribute these 
charge transitions to the so-called corner quantum dots, 
formed in the top corners of the nanowire where the elec¬ 
tric field is maximum, when the channel is in the sub¬ 
threshold regime [T5]. These dots have relatively large 
charging energies, E c ^ 18meV, and a strong top gate 
lever arm, a tg ~ 0.85, similar to what has been reported 
elsewhere mans]. However, one single transition, indi¬ 
cated with a black arrow in Fig. [ljc) , is more coupled 
to the back-gate than the other charge transitions: this 
corresponds to a localized state lying further away from 
the nanowire-top gate interface. 

Such subthreshold resonances are commonly at¬ 
tributed to single donors located in the channel mm, 
which is likely to be the case in the present sample, as 
the charging event occurs very close to the threshold in¬ 
dicating that the impurity state is close to the conduction 
band edge, as expected for phosphorus atoms in silicon 
[21] . In addition, from the stability diagram and the 
back-gate/top-gate lever arm ratio, we can place the im¬ 
purity far below the surface, eliminating any interface 
charge trap or electron puddle, while in the same time 
the number of bulk defects in the nanowire active re¬ 
gion (30 nm x 30 nm x 11 nm) is a few orders of magni¬ 
tude lower than the average number of dopants that were 
implanted (~5 dopants in the channel). Finally, no sub¬ 
threshold signatures were observed in any of the undoped 
devices (see Ref. [15]). 


III. CHARGE DYNAMICS 

A. Interdot charge transition and quantum 
capacitance 

In Fig. |2^a) we focus on a particular region of inter¬ 
est in the charge stability diagram, showing the classic 
signature of a double quantum dot through the presence 
of an extra ridge at the intersection between the charge 
transitions of a corner dot and the donor. The charge 
number (1,1) corresponds to one electron located on the 
corner dot and one at the donor site, and (0,2) indicates 
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FIG. 2. (a) A close-up of the charge stability diagram shows 
the charge transition between a corner dot and a donor, with 
charge occupancies labeled as (corner dot, donor). The * indi¬ 
cates the avoided level crossing sketched in (b) the energy lev¬ 
els of the hybrid double quantum dot, where the quantum dot 
and the P atom form bonding, | g), and antibonding |e) states 
when the detuning, e, is near zero. For simplicity, spin effects 
are omitted here and will be discussed later. For |e| ^>0, the 
double dot has ionic-type wavefunctions where the charges are 
localized on the dot or donor, (c) States | g) and |e) have oppo¬ 
site values of quantum capacitance and hence (d) monitoring 
the interdot charge transition (*), as a function of applied 
microwave frequency, shows a reduction in signal amplitude 
upon resonance with the |gr):|e) transition (see S3 for more 
details). 


that both electrons are on the donor. This charge assign¬ 
ment is deduced from two complementary measurements. 
First, we exploit the phase contrast of the reflected signal 
on the stability diagram (see supplemental information, 
Fig. SI), to assign the two corner dots’ charge occu¬ 
pancy. We deduce that the corner dot transition, seen in 
Fig. [2|a), corresponds to a transition from (0,1) to (1,1). 
Second, from a magnetic field dependence of the donor 
charge transition, we attribute the charge degeneracy to 
a D°:D~ transition (see supplemental information, Fig. 
S2), where D° is the neutral donor and D~ the anion. 
Moreover, we observe this transition very close to the 
threshold, as expected [22] . 

We now investigate the coupling between the donor 
and one of the corner dots. The reflectometry signature 
of a charge tunnelling between the donor and a corner 
dot is obtained by measuring the resonator response at 
an intersection between their charge transitions. The sig¬ 
nal amplitude A at the ICT can be shown to be a func¬ 
tion of the quantum capacitance, its model is described 


in ref. [23] and depicted in Fig. |2|b,c). When the energy 
levels of the donor and the corner dot are brought into 
resonance, the tunnel coupling between the two quan¬ 
tum systems gives rise to a set of molecular orbitals: a 
bonding and an anti-bonding state. The quantum ca¬ 
pacitance is directly proportional to the curvature of the 
eigenenergies with respect to detuning [16i, 

rP 1 E 

C q = -(ea ) 2 ^, (1) 

where a is the coupling between the resonator and the 
double quantum dot, E the eigenenergies and e the de¬ 
tuning energy. The quantum capacitance is then maxi¬ 
mum for e = 0, where the curvature of the eigenenergies 
is maximum. Notably, the quantum capacitances of the 
bonding state, | g), and antibonding state, |e), are of op¬ 
posite sign. 


B. Tunnel coupling and charge dephasing time 

We exploit these quantum capacitance signatures to 
characterize the tunnel coupling between the donor and 
the corner dot as well as the charge dynamics in the sys¬ 
tem, starting by performing MW spectroscopy using a lo¬ 
cal antenna coupled to the top gate. Fig. Id) shows the 
relative change of signal, A/Aq, at the ICT as a function 
of the MW frequency, where A 0 is the signal amplitude 
without any MW applied. The MW excitation causes a 
fraction of the ground state population to be promoted to 
| e), reducing the averaged quantum capacitance. As a re¬ 
sult, the maximum change of signal occurs when the MW 
frequency matches the tunnel splitting, which in our case 
gives A = 2 t c = 5.5 GHz. We use a Gaussian function to 
fit the data, from which we extract the charge dephas¬ 
ing time T£ ~ 200 ps, very similar to the one measured 
in a double-donor system in the same type of nanowire 
transitors [20] . The present dephasing rate is close to 
the tunnel splitting, therefore implementation of quan¬ 
tum protocol that uses coherent charge tranfer to medi¬ 
ate spin information, such as coherent transfer by adi¬ 
abatic passage (CTAP)[2T, would have limited fidelity. 
However, while the tunnel splitting might be tunable in 
alternative device architectures [20] . the dephasing time 
may be improved either by going to lower temperature if 
it is phonon-limited or by improving the charge stability 
of the device if it is charge noise-limited. 


C. Charge relaxation 

To confirm that the dephasing rate is not limited by 
phonon-induced charge relaxation, we present measure¬ 
ment of the charge relaxation time, T\, between the bond¬ 
ing and anti-bonding states. We use MW excitation to 
populate \e) and measure how fast it decays to the ground 
state. The relaxation time is expected to be short com¬ 
pared to our demodulator and voltage amplifier band- 
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FIG. 3. (a) Interdot charge transitions measured under con¬ 
tinuous resonant MW (light blue), chopped MW at 10MHz 
(dark blue) and at 10 kHz (purple), (b) Renormalized av¬ 
erage charge state population difference as a function of the 
chopping period r. The data are fitted using Eq. [2] to give 
the charge relaxation time T\. 


width (^ a few MHz), making a transient measurement 
impossible in our case. Instead, we use a procedure de¬ 
veloped by Petta et al. [25: , where the MW excitation is 
chopped at some frequency 1/r with a 50% duty cycle, 
while time averaging the signal at the ICT. We define the 
charge polarisation as P — (P g — P e ) with P g and P e the 
ground and excited state populations. 

When the period r is long compared to the time 
averaged polarisation is (P) ~ 1/2 (P sa t + Pth)- This is 
because the polarisation is, to a good approximation, at 
saturation, P sa ti during the first part of the cycle (MW 
on) and at thermal equilibrium, P^, during the second 
part (MW off). When r is short compared to Ti, the 
system has no time to relax to the ground state and it 
is then at saturation throughout the cycle, giving (P) « 
P sa t • Charge relaxation takes the system between these 
extremes, giving: 

(P) = — + — + ( P sat ~ Pth) — ~ ^ T/2Tl) ■ (2) 

In the present study, the normalized amplitude A/Aq 
at the ICT is directly proportional to ( P) /P t h , thus 
P sa t/Pth is obtained from the ICT signal under contin¬ 
uous MW excitation (see Fig. |3ja)). The power depen¬ 
dence of the saturation level is plotted in the supple¬ 
mental information, Fig. S3. Also shown in Fig. Up) is 
the ICT under MW excitation chopped with 100 ns and 
100 /is time periods. As expected, the signal amplitude 
is greater when the charge has more time to relax. 

The time averaged polarisation recorded for different 
values of r is presented in Fig. [3jb) and is fitted to give 
Ti ~ 100 ns. Since T\ T 2 , charge coherence is not lim¬ 
ited by coupling with phonons. Therefore, background 
charge fluctuations or noise in the gate voltages are ex¬ 
pected to be the main dephasing sources as they induce 
fluctuations in both tunnel coupling and detuning en¬ 
ergy [26] . 
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FIG. 4. (a) The quantum capacitance signature of the inter¬ 
dot charge transition (ICT) become strongly suppressed with 
increasing magnetic field, as the T_ triplet state becomes the 
ground state, as shown in (b) the energy level diagram near 
the (1,1) to (0,2) transition in the presence of a magnetic held. 
esT denotes the value at which the singlet and T_ states inter¬ 
sect, and is responsible for the shift and asymmetry observed 
in (c), the ICT measured at OmT, 100 mT and 200 mT. 


IV. SPIN BLOCKADE 


We now investigate spin-related effects in the system, 
in order to demonstrate unambiguously that the ICT is 
of even parity. Fig.[4ja) displays the ICT as a function of 
V tg and magnetic held, B. It shows that the reflectometry 
signal disappears with increasing magnetic held and that 
above ^ 1.2 T the variation of quantum capacitance at 
the ICT has completely vanished. This may be explained 
by the fact that i) at e = 0, as the magnetic held is 
increased, the population of the singlet state decreases 
while the population of the T_ triplet state increases 
(and eventually dominates for qiibB ksT); and ii) 
the triplet state has a linear dependence on detuning and 
hence zero quantum capacitance signature according to 
Eq. [T] A similar response has been observed in an InAs 
double quantum dot coupled to a rf resonator [23] and 
has been theoretically investigated in Ref. [27 . 

The shift in the maximum amplitude position observed 
at low held (see Fig. [4j[c)) can then be understood by con¬ 
sidering the detuning value at which the singlet and T_ 
states intersect, esr , which shifts with magnetic held, in 
accordance with Fig. |4^b). For e < cst , there is a singlet 
ground state and the quantum capacitance has a hnite 
value, albeit one which decreases for more negative cst • 
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For e > €st there is a T_ ground state and the quan¬ 
tum capacitance vanishes. The asymmetric lineshape ob¬ 
served is a signature that the system has a significant 
population which follows the ground state, (i.e. at the 
S:T_ intersection) and may be understood by the pres¬ 
ence of an avoided level crossing combined with the fact 
that these measurements are inherently multi-passage ex¬ 
periments. The mixing between these two states is gener¬ 
ally induced by the magnetic field gradient created by the 
nuclear spin bath [9]. In contrast, in this silicon hybrid 
double dot, the mixing should be dominated by the hy- 
perfine interaction with the nuclear spin localized at the 
donor site. As a result, neglecting the Overhauser field 
created by the 29 Si, the electron located on the donor is 
detuned by S = AS.I where A is the hyperfine coupling 
constant, S the electron spin and I the donor nuclear 
spin. In the case of phosphorus, 1=1/2 and A = 117 MHz, 
the S(l,l)-To coherent evolution at very low exchange 
coupling could be driven at ~ 50 MHz. In the case of 
bismuth, 1=9/2 and A = 1.48 GHz, the driving frequency 
should be even higher ranging from 300 MHz up to 3 GHz 
depending on the nuclear spin state. 

V. OUTLOOK AND CONCLUSION 

In conclusion, we have investigated a hybrid double 
dot system formed by a single donor and a corner dot 
in a single silicon nanowire transistor. Combining rf- 
reflectometry with microwave spectroscopy has allowed 
us to determine a tunnel splitting 2 t c = 5.5 GHz as well as 


characterize the charge dynamics. The charge dephasing 
rate 5 GHz is similar to the tunnel splitting, and 

could be reduced by removing the nitride spacers in the 
device, known to possess a large trapped charge density. 
Also the tunnel coupling might be tuneable in a differ¬ 
ent architecture, with a split top-gate for instance, such 
as that used in Ref. [20] . Furthermore, we have demon¬ 
strated spin blockade at the interdot charge transition 
due to the presence of singlet and triplet states. The 
time evolution of such a singlet-triplet qubit should be 
governed by the donor nuclear spin, enabling controlled 
rotation gates when combined with NMR excitation f28j . 
Finally, exploiting the interdot exchange coupling would 
allow a SWAP operation between the corner dot and the 
donor atom spin state and eventually its storage in the 
donor nuclear spin using electro-nuclear double resonance 
techniques. 
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Supplementary information for “Charge dynamics and spin blockade in 

a hybrid double quantum dot in silicon” 


I. STABILITY DIAGRAM CHARGE ASSIGNMENT AND MAGNETIC FIELD DEPENDENCE 


RF reflectometry readout is sensitive to the charge dynamics in the double quantum dot and allows us to map 
out the stability diagram. We can extract the phase response of the device from the I and Q channels. Figure 
[ST] shows the charge stability diagram of the donor-dot system, with the phase response data plotted on the right. 
The signal intensity in a single channel does not vary enough to allow us to differentiate between the systems under 
study. However, the phase response data shows a strong variation in intensity between the different quantum dots in 
the channel. The degree of visibility of its resonances in the phase response should be characteristic of a particular 
quantum dot, since it is related to the proportion of its signal in the dissipative versus dispersive regimes. Using this 
method of discrimination, we are able to isolate three sets of lines in Figure SI (a); a corner dot uncoupled to the 
donor (blue dotted line), a corner dot that couples to the donor (red dotted line) as well as the donor charge transition 
(green dotted line). 



FIG. SI. (a) Numerical derivative, filtered, of the amplitude response as a function of V t g and Vb g - The last digit of the charge 
assignment corresponds to the donor atom (green dashed line), the second digit to the coupled corner dot (red dashed line) 
and the first digit to the uncoupled corner dot (blue dashed line). The unassigned transitions are attributed to another donor 
in the channel, (b) Phase response of the coupled dot-donor system as a function of Vtg and Vb g - The visibility of the phase 
response allows discrimination between the two corner dots. 


Figure S2 shows the magnetic field dependence of the interdot charge transition (ICT) at e = 0 as well as of the 
dopant line. The interdot quantum capacitance signal has a strong field dependence, explored in the main text. As 
observed in other studies [23|, the quantum capacitance signal under Pauli Blockade rapidly diminishes in intensity at 
this even parity transition under increasing magnetic field, and it also shifts to more negative values of detuning. This 
allows us to assign the (1,1) and (0,2) states based on the shift in detuning since the T_-singlet crossing shifts in the 
direction of the (0,2) state. Additionally, while the charge assignment in the stability diagram commences from zero 
electrons for clarity (at a point in top gate voltage below which we do not see any more dot transitions), in essence 
we can assert only that we are at an even parity transition at the ICT. 


The dopant line does not diminish under field, but it exhibits a gradual shift towards higher top gate voltages. In 
general from this we can infer the spin polarity of the tunneling electrons: a shift towards higher top gate voltage, or 
higher energy, indicates tunneling of spin-up electrons. In our system, as studied previously EUHS], this shift towards 
higher energy is consistent with a D - state, in which electrons form a two-electron singlet. Finally, the D + :D° charge 
transition of the donor is not visible in our system. This may be explained by a small tunnel rate between the donor 
and the reservoir. Indeed, the D + :D° is deeper below the conduction band than the D°:D“, and as a result, the 
tunnel barrier between the reservoir and the donor is more opaque [29] . Due to the fact that our RF reflectometry 
technique is sensitive to charges tunneling at a frequency comparable to the excitation [30] , we are unlikely to see the 
D + :D° charge transition. 
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FIG. S2. (a) Resonator response as function of of Vt g and magnetic field B , and (b) as of of Vtg and Vb g - The two charge 
transitions of interest are labelled with stars in (b). 


II. MICROWAVE SPECTROSCOPY POWER DEPENDENCE 


Figure [S3] shows the interdot charge transition under a continuous 
power and detuning. 


microwave excitation of 5GHz as a function of 
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FIG. S3. (a) Resonator response as a function of \ tg and the relative microwave power broadcasted by the antenna, (b) 

Amplitude of the interdot charge transition as a function of the microwave amplitude. Equation 0 is used to fit the data. 

In a two level system, Landau-Zener-Stiickelberg (LZS) interferometry theory provides us with the respective slow 
passage (Av < A 2 ) and fast passage (Au A 2 ) conditions, where A is the driving amplitude, v is the frequency, and 
A is the energy splitting between the excited and ground states. As we vary the power we supply to the microwave 
antenna, we are varying the driving amplitude A , which in turn takes us from the slow to fast passage regime. This, 
as shown in Figure |S3[ gives rise to the expected squared Bessel function behaviour at zero detuning for multiple 
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passage, as theoretically investigated in Ref. m- The MW driving tunes the charge polarisation as: 

P g -P e = J 2 0 (a) ( 1 ) 

where Jg(<a) is the square of the 0-th order Bessel function evaluated at a = eV ac /hv. We send microwaves through 
an antenna close to the sample, which is then assumed to vary our gate voltage by V ac - Our data, while fitting well to 
the Bessel function, also allows us to calibrate the power arriving at the sample by measuring the distance between 
the zeros of Jg, occurring at aq = 2.4 and ot 2 = 5.5. Figure S3(b) shows the calibrated result. 

While we do see a characteristic power dependence at detuning e = 0, no higher order photon assisted tunneling 
has been observed in the high frequency limit ( hf > 2 t c ). This can be explained by considering spin effects in the 
double dot at the (0,2) to (1,1) transition: for negative detuning, the ground state is the singlet S(0,2) while the 
excited state is a mixture of triplet states (To, T_ and T+) and singlet S(l,l). As a result, a MW induced transition 
between these two states requires a spin-flip process (Am = 1). Hence, in the absence of spin orbit coupling [32] . or 
a particular electron-phonon interaction [33] . MW photon absorption is forbidden. 


III. SUPPLEMENTARY DEVICE 


Figures S4[a) and (b) present measurements done on a similar device as the one presented in the main text, but 
with a wider and larger channel (240 nm x 100 nm x 10 nm). A tank circuit with a resonant frequency of 671 MHz is 
connected to the top gate. A small AC excitation is applied to the drain at frequency of 1.1 kHz. The amplitude of 
the resonant side band (671.0011 MHz) is measured as a function of top gate and back gate using a spectrum analyzer. 
The presence of numerous charge transitions in Fig S4'a) is attributed to the larger number of dopant in the channel 
compared to the device presented in the main text. Fig. [S4^b) shows a close-up of the stability diagram where a 
donor-based interdot charge transition is visible, similar to what is observed in Fig 2(a). 



-2 0 2 0.50 0.57 0.64 


V 


bg 


(V) 



(V) 


FIG. S4. (a) Modulus of the derivated reflected amplitude as a function of Vt g and Vbg- The white rectangle corresponds to 

(b) a close-up of the stability diagram showing a strong signal based on a interdot charge transition. The extra line crossing 
the interdot charge transition is attributed to an uncoupled corner dot. 











